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Abstract 

The Born approximation to the quark-gluon-exchange mechanism for K+N scat- 
tering is used as a starting point to generate a potential for this system. The valence 
quark wave function of the nucleon is generalized from a single Gaussian to a sum 
of Gaussians in order to have a more flexible representation than previous work. 
We obtain a potential derived from a valence density given by lattice calculations. 
By comparing with a recent amplitude analysis it is found that the strength of the 
quark-gluon based potential needs to be increased by a factor of order 2-4 relative 
to the normalization given by more traditional values of the governing parameters. 
The method is used to estimate the change in effective K"''N amplitudes which would 
result from changes in the valence quark distributions or strength of the interaction 
which might arise from nuclear medium effects in scattering from nuclei. 

PACS numbers: 13.75.Jz,25.80.Nv,12.39.Pn 
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I. INTRODUCTION 

The K"*" meson is useful in the study of the strong interaction because of its simple struc- 
ture, i.e. it has no light valence antiquark. This feature means that no 3-quark intermediate 
state can exist, thus removing the most common resonances from the interaction. Although 
the interaction is weaker than other typical strong interactions it is still substantial and must 
be treated as such. 

One of the principal motivations of this work is to obtain the relationship between the 
distribution of the constituents of the nucleon and the K^-N scattering amplitudes. In Ref. 
[l| it was suggested that the increase in size of the nucleon in the nucleus could be observed 
as an alteration in the scattering of the kaons from nuclei. In that case the model used to 
connect the change in size with scattering was a strongly repulsive well whose radius was 
allowed to vary. While this simple model may be expected to encompass much of the basic 
physics, clearly a model based on fundamental degrees of freedom is to be preferred. 

Barnes and Swanson [2!] treated the K^-nucleon interaction due to the exchange of one 
gluon and two quarks, one of the lowest order process dealing directly with the underlying 
degrees carried out a similar program for meson- 

meson scattering but we do not address that sector here). Their calculation of a lowest 
order exchange must be iterated to obtain a full evaluation of the scattering amplitude. A 
common way to do this is to treat their result as the Born approximation corresponding to 
a potential and use a wave equation to generate an approximation to the full amplitude. 

Such a program has been very successful in the case of one pion exchange. The one-pion- 
exchange potential is able to predict many properties of the deuteron 0, [sl, . In the pion 
case the exchange of a single particle gives rise to a local potential. In the present (quark- 
gluon exchange) case the potential has one term which is local and three others which are 
non-local. 

In the case of one pion exchange, a potential is essential in order to be able to produce 
the deuteron bound state. In the present work the potential is needed to calculate realistic 
phase shifts for comparison with data. This potential will then enable us to relate possible 
changes in valence quark distribution or strengths to the changes in phase shifts which would 
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result. 

Because they included no mechanism that would lead to spin dependence, the results of 
Ref. 2j] were necessarily spin-independent. The experimentally determined partial waves 
with ^ greater than zero show a strong spin dependence [3] so that a comparison with these 
partial waves is not possible and they concentrated their work on the £ = partial wave. 
We will do the same. 

n 

In a later, similar, program N. Black|8l] calculated the higher partial waves with quark- 
gluon exchange including the spin dependence. For T=l he found qualitative agreement 
except for the P3/2 partial wave where the sign was wrong. For T = the phases shifts 
were considerably smaller than the experimental ones, although there was a tendency for an 
alternation of sign similar to what is observed. In any case, the £ = partial wave dominates 
the isospin one amplitude at low to moderate energies. 

The Barnes and Swanson{^ calculation assumes a Gaussian wave function for the distri- 
bution of the quarks in the nucleon. This is a rather restrictive assumption on the shape 
and does not represent the functional form found in lattice calculations 9|, llO| . To repeat 
their calculation directly with more realistic wave functions would be difficult because of the 
9-dimensional integrals involved. With Gaussian wave functions these integrals can be done 
analytically so that an alternative method of treating this problem is by representing the 
desired wave function as a sum of Gaussian functions. In this case the integrals can be done 
term by term but involve cross terms in which the Gaussian parameters in the initial and 
inal states are different. Such an approach requires a derivation of the expressions of Ref. 
2! for this more general condition. We carry out this step and then calculate a potential 
for a wave function consisting of the sum of two Gaussian functions. One could include 
more terms but at this stage the knowledge of the desired wave function is probably not 
sufficiently accurate to warrant the effort. ^ 

In the calculation of Barnes and Swanson^] the Born approximation to the amplitude for 
K+N scattering is given by 
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where 

F.(k,kO=e--'=^+^'^'^', mn = ^^, (2) 

k and k' are center of mass momenta, = |k| = |k'|, 7 = as/m^q-, En is the total energy of 
the nucleon and Ek is the total energy of the kaon in the center of mass, is the running 
strong coupling constant and ruq is the light (constituent) quark mass. Expressions for the 
Oj, 6j, and rji are given in Ref. [3] for simple harmonic oscillator wave functions. They 
are functions of p, the ratio of light to strange quark (taken to be 0.6), a, the parameter 
governing the size of the nucleon and (3, the parameter determining the size of the kaon 
through the ratio g = 

They considered four diagrams for quark and gluon exchange. For the first diagram the 
momentum transfer is given by the two exchanged quarks with gluon exchange between 
them so the resulting interaction is local as indicated by the fact that ai = 61. For the other 
three diagrams this is not the case and the interaction is non-local. Barnes and Swanson 
calculated an equivalent local potential valid at zero energy only. 

The isospin dependence is contained in the weights, Wi, which are given by 

W={o,i.0.1}:r = Oa,>d M = {l.l.i.l}:r=l. (3) 

One sees that there is a bias toward the isospin unity amplitude being larger than isospin 
zero. If the contributions from each of the four diagrams were equal (they need not be, of 
course, and are not) then the T = 1 amplitude would be 7/3 of the T = amplitude. 

It is interesting to examine the low-energy limit of this amplitude. In this case the 
functions, Fj(k, k'), become unity. The functions, rji depend only on g so, in this limit, the 
absolute size of the kaon and nucleon systems does not enter but their relative size remains 
very important. 

In Section II we derive expressions for the potential for arbitrary values of the parameters 
in each term, in Section III we study the results with the use of various single Gaussian 
wave functions for the valence quarks, and in Section IV we develop the two-Gaussian form 
for the wave function. In Section V we give a summary of the formulas needed to calculate 
the parameters for the two-(or multi-) Gaussian potential, in Section VI we show the basic 
fit with the two-Gaussian form and study the consequences of a variation of the radius or 
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strength (as for the case of immersion in nuclear matter). Section VII gives the results of 
a calculation of the off-shell amplitude which results from the two-Gaussian form for the 
potential and Section VIII gives some conclusions. 



II. DEFINITION OF THE POTENTIAL 



The potential will be given in terms of the amplitude by 



nr,r') = — ^^^.^.F.(r,rO (4) 
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Fi{r,r') = — / dkdk'e^^''-'^''''Fi{k,k'). (5) 



where 

To get a fully off-shell amplitude which reduces to this form on shell but is a function 
only of momentum transfer if a = 6, we write (dropping for the moment the index i) 

F(k,k') = e-t(^^+^")+i'-'^''' (6) 

where k and k' are independent and are assumed not to be constrained by the on-shell 
condition. We hasten to point out that this is a choice and that other dependences are 
possible. We make no claim to uniqueness. We now choose new vector variables as 

= ^ ^ ^ ; q = k — k' so that \s.— k + ^; k' — K — ^ (7) 

allowing us to express the form as 

For the special case of a = 6 we can take the Fourier transform on q to get the local 
potential 

3 

Fi(r) = y (iqe^'i '"e-^»i«' = (^^^ ' e-^'/«i. (9) 
In the general case we can write the Fourier transform as 

F(r,r') = ^ y dk(ik'e-^'"'-e^'''-'-'e-i(^'+^")+''-'^''' (10) 
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(27r)6 

We use the formula 

/ dre'^^ 'e-'"'''^ = '—e'i^ (12) 
J a^^ 

to do the two Gaussian integrals separately to obtain 



^ 2 



F(r, r') = ^ V-^e (13) 

(27r) - 



A , 2 , /2, CXD 

47r -Mi+n) 



e 



2{a^-b^) 



(27r)3v/^^352-^ ^ 



5^yr(^)>^r(nQ.(^) (14) 



where 



g,(^) = {-iYji{iz) (15) 



and je{iz) is the spherical Bessel function of imaginary argument. The local potential can 
be recovered by taking the limit 6 — a in the non-local expression. 
The wave equation to be solved is 

V'i^iv) + ^I4(r)V^(r) + dv%{v,T')i,{v') = kli,{v) (16) 

j=2,3,4 

from which the solution for the wave function and phase shifts can be obtained by standard 
matrix techniques. 

Because local potentials are relative easy to deal with, a commonly used approximation 
has been to create a local potential which is energy dependent although it is inconvenient 
due to the loss of orthogonality of wave functions at different energies. We considered such 
a procedure in this case. Thus we could write 



F(k, k') = e-^(^ +^ )-|('^-'^') , g-(a-6)fco^-|(k-k')^ (^7) 

where fco is the free-space wave number. However, for the standard parameters given in Ref. 







2j for the fourth diagram, the value of 64 is negative so one cannot create a local potential for 
this term in this way. We then abandoned this procedure and continued with the non-local 
potential for the terms beyond the first. Of course, an energy-dependent local potential can 
be generated for the sum of the four terms after the local/non-local system has been solved. 
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a. 


R{N) 




R{K) 


F 


Born 


Potential 


Set Description 


(GeV/c) 


(fm) 


(GeV/c) 


(fm) 




ao,ai (fm) 


ao,ai (fm) 


a Ref. set 


0.40 


0.85 


0.35 


0.49 


1.00 


-0.12, -0.35 


-0.10, -0.23 


b Fit 


0.68 


0.50 


0.43 


0.40 


0.98 


-0.15, -0.31 


-0.11, -0.18 


a' Mod. Ref. 


0.40 


0.85 


0.50 


0.34 


2.06 


-0.16, -0.78 


-0.12, -0.34 


b' Mod. Fit 


0.68 


0.50 


1.15 


0.15 


4.84 


-0.21, -1.89 


-0.11, -0.31 


c Lattice 


0.61 


0.56 


0.93 


0.18 


3.36 


-0.18, -1.32 


-0.11, -0.31 



TABLE I: Parameter values for the single- Gaussian fit to the quark density. R{N) is the rms value 
of the inter-quark distance in the nucleon and R{K) is the sanie quantity for the kaon. Also given 
is the original reference set as well as the fitted set of Ref. [2^. Sets a' and h' are the result of 
holding the nucleon size fixed at the unprimed values and adjusting j3 and F to fit the phase shifts. 

III. RESULTS FOR SINGLE GAUSSIAN 



Many different assumptions have been made about the spacial extent of the valence 
quark distribution of the nucleon so that we need to consider a possible range of models. 
One approach which should be useful is the calculation on the lattice. We will use that 

n n 

of Lissia et al. [9[ although the results of Alexandrou et al. [10] are similar. Figure |^ shows 
(dotted line) the result of a fit of single Gaussian to the data from Lissia et a/. [9]. For 
this single-Gaussian case we can use the equations from Ref. [s] directly to compute the 
amplitude. Since the light quark mass and the strong coupling constant enter only in the 
combination, 7 = a^/m^, rather than consider changing these quantities independently we 
calculate relative to a standard value of 70 = 0.6/0.33^ = 5.51 (GeV/c^)~^. We define 



F = 7/70 anc 
Reference 



<:eep the ratio of light to strange quark fixed at p = 0.6. 
2| first calculated with a "reference" set of parameters taken from quark sys- 
tematics. For this "reference" set we find that, although the scattering lengths for the Born 
approximation are in ag , the potential solution is not. 

The phase shifts from this calculation are shown in Fig. O Although the Born result has the 
correct ratio of iso-triplet to iso-singlet scattering length, the potential result does not. Since 
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FIG. 1: One and two Gaussian fit to the proton wave function. The points are from Lissia et oZ.[9|]. 



this ratio is sensitive to the relative hadron sizes we can change the kaon size parameter to 
f3 = 0.5 GeV/c to get the correct ratio in the potential calculation and increase the strength 
of the potential by a factor of 2.06 to get the right magnitude. In this manner a good fit for 
the potential calculation is found (see case a' in Table [T]). 

If we use the "fitted" values obtained by Ref. ^ we find case b in Table [H For this value 
of the nucleon size we need to alter the strength of the potential by a factor of 4.84 and 
choose (3 = 1.15 GeV/c to get the correct values for the scattering lengths (Table H]). 

If we fix the value of a = 0.61 GeV/c from the single Gaussian fit to Lissia et ai, increase 
the strength of the potential by a factor of 3.36 and take [3 = 0.93 GeV/c we are able to 
fit the scattering phase shifts well (case c in Table [T]). Figure [3] shows the result of this 
calculation along with the "reference" set of Ref. Q| . 

Figure H] show the results for the individual diagram potentials without isospin factors 
from Eq. [3] (equivalent local potential at zero energy for all but the first) for the "reference", 
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FIG. 2: Results of Barnes and Swanson with their reference and fitted set in Born approximation 
(dashed Unes) and the corresponding phase shifts from the solution with the potential (solid lines). 
The phase shifts from Ref. 71] are shown by the solid dots. 

"fit" and lattice input of case c. Also shown (short dash curve) is the result of changing a 
from 0.61 GeV/c to 0.55 GeV/c {Rrms = 0.56 fm to Rrms = 0.62 fm). 

It is seen that the ranges of the potentials from the second and fourth diagrams (the only 
ones which contribute to the T = amplitude) appear to decrease with increasing nucleon 
radius whereas the range of the potential from the first and third (which dominate the isospin 
one amplitude) increase. One can trace this behavior back to the quantities rji, ai, and 6j. In 
fact, the "true" range parameters (a^, 6j) always increase with increasing Rrms but ri2 and 
7]^ decrease rapidly. 
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FIG. 3: S-wave phase shifts using a wave function with a single value of a obtained from the fit 
to Lissia et al. [91]. Born approximation (short dashed lines) and the corresponding phase shifts 
from the solution with the potential (solid lines). The result of the reference calculation and its 
potential corrected values are shown by the long dash and dash dot curves respectively. The phase 
shifts from Ref. 71] are shown by the solid dots. 

Figure [5] shows the full potentials for the two isospins for case c. It is seen that the 
potential is very near to a hard-core potential for case c. This result can no doubt be 
traced to the fact that for both isospins the phase shifts are nearly linear as a function 
of the center-of-mass momentum as is the case for a hard-sphere interaction. Also shown 
(short dash curve) is the result for a = 0.55 GeV/c. As expected from the remarks above, 
the isospin zero potential has the opposite variation with nucleon size as the isospin one 
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FIG. 4: Potentials without the isospin factors for K^N scattering for three of the parameter sets 
given in Table [H The upper left hand panel shows the local potential corresponding to the first 
diagram and the rest show the corresponding equivalent local potentials at zero energy. The short- 
dash curve with c shows the effect of changing a from 0.61 to 0.55 GeV/c. 

potential. 

These results show that the Born approximation (especially for T = 1) is not accurate. 
The use of a potential (or other unitarizing agent) is essential. The range of the potential 
(related to the size of the hadronic systems) is important. For a more extended potential 
the Born approximation is more nearly correct. 
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r (fm) 

FIG. 5: The full isospin potentials for a single Gaussian at zero energy using the equivalent local 
forms for the non-local potentials calculated with the parameter set c in Table |T] (solid curve) and 
with a = 0.55 GeV/c (dashed curve). 



IV. TWO GAUSSIAN FORM OF THE WAVE FUNCTION 



With the extention of the formulas of Ref. given in the next section one could carry 
out the calculation of a potential using a sum of any number of Gaussian functions for the 
quark wave function of the nucleon. A similar program was carried out for meson-meson 
scattering by Hilbert et a/. [11]. Of course, such a sum would introduce a large number 
of parameters which have to be determined. Since we intend to fit them to lattice data, 
and those calculations are of a limited accuracy at present, we use a sum of only two such 
functions. The wave functions below are expressed in Jacobian coordinates as discussed in 
Appendix A. 
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We take the sum of two Gaussian wave functions, each normahzed independently. With 

^P{p, q) = NaMP^ Q) + NbMP, Q) (18) 
we have for the normalization integral 

j dpdqtP^ip, q) = 21{Nl + N^) + 2NaN, j dpd<iMP, Q)Mp, Q) = 1- (19) 

Since 



! dpdciMP,q)MP,q) = ^X^r f P'dpq'dqe'^' V^-^iJe"^ i^'"^) = 27 



(20) 



the normahzation condition becomes 



(«^ + 

If we take Nb = XNa then 



A^-+ ^°y't' +iV,' = l. (21) 



a 



K + air ^ 



1. (22) 



We are free to take X, aa and ab as we wish and, as long as this condition is satisfied, the 
wave function will be normalized. 



V. SUMMARY OF FORMULAS 

In order to carry out the calculation of the potential (or the amplitude) for a sum of 
Gaussian forms for the wave function we must evaluate the integrals given in Ref. (Eqs. 
41, 42, 43 and 44 which are our equations IB5[ IB161 IB28I and IB39P for the case where kaons 
A and C have different governing Gaussian parameters (/3) and nucleons B and D also have 
different parameters (a). The evaluation of the integrals is done in appendix B and here we 
give a summary of the results. Auxiliary definitions are given at the end. 



1 1 c| + 

02 = — + TTT^TTT + 



2(l + p)2^2 1252 4(q,2 + 2/3^2 + 

1 1 C2C' 



2(l + p)2^2 12^2 2(a? + 2/3^2 + 

^ ^ [(2 + 7f)(l + p)-672l , ^ [-7r(l + p) + 6(7f + 

3(1 + p) ' 3(1 + p) 

^2 



3a2 8(2a22 + 3/32)' ^ 3^2 4(2^22 + 3/32) 
2[27f + (272"-3)p] . 2 1 

— 3(rT^^^ — ' ^3 = -3(272+1) 



^a6(2a2^ + 3/32)1 
(l-p)2 , 1 , C| + Cf 
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64 = 



4(« + 2/3)(l + p)2 12a2 4(^ + iay 
(l-p)2 , 1 C4CI 



4(d + 2/3)(l + p)2 12a2 2(^ + iay 



2«;2 - t;2 + ^ + (^;2 + ^ ^ + 3^2 _ (^2 ^ ^2)^a 



3(^;2 + «;2) ' 3{v^ + w'') 

q q ^ 2 / 8 \^ 1 



774-p«^«^/3^/3| 



U2+«;2 3 "12 J 
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FIG. 6: S-wave phase shifts with a^, and x chosen from Lissia et al. [9;] and the overall factor 
on the potential chosen to fit the T=l phase shift. The dashed curve shows the result if the value 
of Oa is altered from its fitted value corresponding to the changes in density shown in Fig. [71 The 
phase shifts from Ref. [tI are shown by the solid dots. 



VI. RESULTS FOR A DOUBLE GAUSSIAN WAVE FUNCTION 



Figure [T] shows a fit (solid line) of the two-Gaussian type to the lattice calculation of 
Lissia et al.^\. The values of the double Gaussian parameters are = 0.58 GeV/c, ai, = 
1.32 GeV/c and X = 0.1. This selection of parameters gives a reasonable representation of 
the s-wave phase shift (see Fig. |6]) but a factor of 4.1 has been multiplied by the potential 
to obtain the fit. 
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FIG. 7: The solid curve gives the result for r^/j(r) with Ua, Ofe and X chosen from Lissia et al. 
The curves show the densities that result for various values of with = 1.32 GeV/c and X 
fixed at 0.1. 

The fit with the double Gaussian is shown in Fig. [6] as the sohd fine. Also shown are a 
series of variations of ai (holding a2 and X fixed) used to vary the radius of the system. The 
effect of this variation on the density is shown in Fig. [71 This variation would be supposed 
to come from the partial deconfinement of the quarks in a nucleus. As expected from the 
remarks on the behavior of the potentials as a function of nuclear size for the single Gaussian 
case, the T = phase shift decreases with increasing nucleon size. 

A.*„ugK *e ...o of to.a. cos. sec«o. of nude, to of deute..„ —U Q 
clearly shows a deviation from the predicted ratio with free kaon-nucleon interactions, the 
energy dependence is not in agreement with the simple change in radius estimated in Ref. 
[l| being larger at higher incident momenta. Since the model used to calculate the ratio is 
rather simple, one should consider the possibility that a more realistic model could give the 
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energy dependence necessary. For this reason we compare the change in cross section that 
results from a given change in nucleon radius for two different energies. We also make the 
same type of comparison for changes in interaction strength since that possibility has been 
suggested as an alternative to an increase in nucleon sizejl4|. 

In a series of papers, Caillon and Labarsouque ISj studied this problem with different 
possibilities for the interaction of the kaon with the nucleus. They are able to partially 
resolve this problem but have difficulty predicting the correct dependence on atomic number. 

In Tables HIl llllj iTVl and |V] is shown the effect on the elastic and total cross sections of 
the variation in radius and strength. For each isospin the amplitude 

p2i<5T _ 1 

/. = ^— (37) 

is calculated and the isospin weighting is done. Since at finite energies a linear expansion of 
the exponential does not provide a good approximation and since the variations in the two 
isospin phase shifts with radius and strength go in opposite directions there is not a simple 
linear proportionality. 

From Table [Tll we see that at a beam momentum of about 0.6 GeV/c a change in the 
radius of 15% corresponds to a change in the differential cross section of 14% and a change in 
the total cross section of 18%. A change in the radius of 25% gives corresponding changes of 
25% and 32%. From Table UTTl we see that at a higher beam momentum the same fractional 
changes in radii lead to smaller fractional changes in the cross sections. The trend is opposite 

nn 

to what would be needed to give a better fit to the data[12, ll3j- 

Tables IIVI and |V] show a similar effect for the variation of the cross sections as a function 
of the change in the multiplying factor. It is seen that a significant change (in percentage) of 
the multiplying factor is needed to produce a modest change in the total cross section ratio. 



VII. OFF-SHELL AMPLITUDE 



From the full solution for the wave function one can calculate the scattering amplitude 
for off-shell momenta. We have done so at (essentially) zero energy. Since the value of the 
off-shell amplitude at zero momenta will be equal to the scattering length, dividing this 
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R(fm) 








/o 


fl 


f/i + i/o 


k^a 




ctt/ctt 


0.566 


1.00 


-0.210 


-0.562 


-0.204 + i0.043 


-0.451 + zO.284 


-0.389 + iO.224 


0.202 


1.00 


1.00 


0.652 


1.15 


-0.192 


-0.623 


-0.188 + ^0.037 


-0.474 + iOMO 


-0.402 + iO.264 


0.232 


1.15 


1.18 


0.709 


1.25 


-0.178 


-0.669 


-0.175 + i0.032 


-0.486 + iO.384 


-0.409 + iO.296 


0.255 


1.26 


1.32 


0.744 


1.31 


-0.166 


-0.700 


-0.166 + i0.029 


-0.493 + i0.415 


-0.411 + i0.318 


0.270 


1.34 


1.41 



TABLE II: Dependence of the s-wave amplitude on the quark radius at Pc.m. = 0.35 GeV/c 
(PLab=0.599 GeV/c). R is the rms value of the inter-quark distance in the nucleon. The quantity a 
is the (isoscalar) differential cross section which would come from the s-wave alone and the quantity 
ar is the (isoscalar) total cross section from the s-wave alone. The superscript zero refers to the 
first line in the table i.e. the result obtained from the fit to the data in free space. 



R(fm) 


R/R^ 


So 


Si 


/o 


fl 


|/i + i/o 


k^a 






0.566 


1.00 


-0.312 


-0.816 


-0.292 + zO.094 


-0.499 + i0.531 


-0.447 + iO.422 


0.378 


1.00 


1.00 


0.652 


1.15 


-0.287 


-0.900 


-0.271 + Z0.080 


-0.487 + i0.614 


-0.433 + iO.m 


0.418 


1.11 


1.14 


0.709 


1.25 


-0.267 


-0.963 


-0.254 + i0.069 


-0.469 + iO.673 


-0.415 + «0. 522 


0.445 


1.18 


1.24 


0.744 


1.31 


-0.253 


-1.004 


-0.242 + i0.063 


-0.453 + i0.711 


-0.400 + iO.549 


0.462 


1.22 


1.30 



TABLE III: Dependence of the s-wave amplitude on the radius at Pc.m. = 0.5 GeV/c (Pia6=0.941 
GeV/c). R is the rms value of the inter-quark distance in the nucleon. The superscript zero refers 
to the first line in the table i.e. the result obtained from the fit to the data in free space. 

amplitude by the scattering length one obtains a function which is unity at g = 0. Because 
the wave function is constant at zero energy the result is equivalent to the Fourier transform 
of the sum of the four potentials. 

Figure M shows the result of such a calculation. We see that the half value of the T=l 
function occurs at around 1 GeV/c. T=l dominates nuclear scattering because of the average 
weight in the isoscalar average and the relative smallness of the T=0 scattering length so its 
range is perhaps the appropriate one to use. 
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F 




5o 




/o 


/i 


|/i + 3/0 


cr 


(7/(7° 




4.1 


1.00 


-0.210 


-0.562 


-0.204 + i0.043 


-0.451 + zO.284 


-0.389 + iO.224 


0.202 


1.00 


1.00 


5.0 


1.22 


-0.230 


-0.591 


-0.222 + i0.052 


-0.463 + i0.311 


-0.403 + iO.246 


0.223 


1.10 


1.10 


6.0 


1.46 


-0.248 


-0.618 


-0.238 + i0.060 


-0.472 + iO.336 


-0.414 + i0.267 


0.242 


1.20 


1.19 


7.0 


1.71 


-0.263 


-0.640 


-0.251 + zO.068 


-0.479 + zO.357 


-0.422 + zO.284 


0.259 


1.28 


1.27 


8.0 


1.95 


-0.276 


-0.659 


-0.262 + zO.074 


-0.484 + zO.375 


-0.429 + zO.300 


0.273 


1.35 


1.34 



TABLE IV: Dependence of the s-wave amplitude on the multiplying factor at Pcm. = 0-35 GeV/c 
(PLa6=0-599 GeV/c). The superscript zero refers to the first line in the table i.e. the result obtained 
from the fit to the data in free space. 



F 


F/F° 


So 


Si 


/o 


/i 


3/1 + 3/0 


k^a 






4.1 


1.00 


-0.312 


-0.816 


-0.292 + iO.OM 


-0.499 + zO.531 


-0.447 + zO.422 


0.378 


1.00 


1.00 


5.0 


1.22 


-0.341 


-0.858 


-0.315 + ^0.112 


-0.495 + zO.572 


-0.450 + iO.457 


0.411 


1.08 


1.08 


6.0 


1.46 


-0.367 


-0.896 


-0.335 + i0.129 


-0.488 + i0.610 


-0.450 + io.m 


0.442 


1.17 


1.16 


7.0 


1.71 


-0.388 


-0.927 


-0.350 + i0.143 


-0.480 + i0.640 


-0.447 + i0.516 


0.466 


1.23 


1.22 


8.0 


1.95 


-O.iOG 


-0.955 


-0.3G3 + /0.15G 


-0. 172 + /O.GGG 


-0.111 + /0.538 


0.187 


1.29 


1.27 



TABLE V: Dependence of the s-wave amplitude on the multiplying factor at Pcm. = 0.50 GeV/c 
(PLa6=0.941 GeV/c). The superscript zero refers to the first line in the table i.e. the result obtained 
from the fit to the data in free space. 

VIII. DISCUSSION AND CONCLUSIONS 

We have calculated the kaon-nucleon scattering amplitude as a function of the nucleon size 
and interaction strength in the one-gluon-exchange approximation. Two interesting features 

are observed: 

1) The isospin one and isospin zero amplitudes do not behave in the same manner under 
a change in nucleon size. The isospin zero amplitude actually decreases with increasing 
nucleon size. Although the potential used is not unique, the fact that the behavior is due 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

q (GcV/c) 

FIG. 8: Ratio of the ofF-shell to on-shell amplitude at zero energy for T=l (solid) and T=0 (dashed). 
The function falls to 1/2 for T=l at about 1 GeV/c. 

to the sharp decrease of 172 and 774 with increasing nucleon size leads us to beUeve that the 
behavior is more general 

2) A significant renormahzation of the potential is required to fit the experimental am- 
plitudes. The factor of 4.1 on the strength (relative to a fixed value of 7 = as/m^q) of 
the interaction necessary to bring the quark-gluon exchange into agreement with the s-wave 
phase shifts is considerable. We have found that the factor needed depends on the extent 
of the quark distribution of the nucleon and kaon. If this result is confirmed then we must 
expect additional quark diagrams to play an important role or, alternatively, the exchange of 
color singlets (for example two-pion exchange with excitation of the K*) may be important. 
The very small range of the kaon distribution is also disturbing. Because the potential which 
results from the fit is very near to that of a hard-sphere, one might think that Pauli blocking 
might play an important role ^\ ■ 
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APPENDIX A: JACOBY COORDINATES 

The integrals needed for the calculation are more easily done if the momentum variables 
are expressed in terms of Jacoby coordinates. Barnes and Swanson use a nucleon wave 
function represented as 



, 3 



34 _ P1+P2+P3-P1-P2-P2 P3-P3 P1 

(Pi,P2,P3) = ^-e ^ . (Al) 



Expressed in terms of coordinates 



P = Pi-P2; Q^^i^; q = P3-Q (A2) 



the wave function becomes 



Kj ^ P y 

p, q) = — 3 — e~4^e~3^. (A3) 



We can transform these wave functions into coordinate space. In the Fourier transform 
we have 

Pi ■ ri + p2 ■ rs + P3 ■ rg = Q ■ R + ■ r + q • rg (A4) 

where 

R = n + r2 + ra; r = ri - r2. (A5) 
Thus the dependence in r-space is 

0(r, rs) oc e 4 g 4 (^A6) 



givmg 

< >= R\N) = (A7) 

for the expectation value of the square of the inter-quark separation in the density. 
The normalization in this form follows from the normalization given {2] 

J C?PiC/p2C/p30^(pi, P2, P3)5(Pc.m. - Pi " P2 " P3)<5(Pc.m. " Pi " P2 - Ps) (A8) 
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= I dpdQrfq02(p, Q)^(P,^ -3Q-q)^(P;.^ -3Q-q) = ^^^""^'^^ ^^'"■^ J dpdci<j>\p, q). 
The integral over p and q should give 27. 

3 

j dpdq(f)'^{p,q) = {^T^f^-^j p^dpe~^ j q^dqe'^ = 27. (A9) 

The kaon wave function is written in terms of the relative momentum as 

1 _ 

<Pkaon{p) = ^-re ^. (AlO) 

APPENDIX B: DERIVATION OF THE AMPLITUDE EXPRESSIONS 

We calculate the expression of Ref. |2] for the case where the initial and final Gaussian 
parameters are not the same in order to enable a representation of the wave functions as a 
sum of Gaussian functions. 

The strategy is the following: 

a) Rewrite the nucleon wave functions in Jacoby coordinates in the form of Eq. IA3[ 

b) Transform to r-space to get the proper combinations of the momenta. 

c) Perform the momentum integrals to find 3 (5-functions in r 

d) Carry out the remaining 3 Fourier transforms. 

One has two quantities to follow as the integrals are done: 

1) The imaginary part of the exponent and 

2) The real parts of the exponent which represent the Gaussian functions and, when the 
integral is completed, give the dependence of the amplitude in the form of a,, and 6j. 

In this process the factors in front of the exponentials are ignored. They are obtained 
most easily by setting k = k' = and performing the momentum space integrals directly. 
The results for these normalizations {rji) are given in the summary. 

The basic Fourier transform used is given by Eq. [121 A useful transformation is 

obtained with the change of variables 

t = r + ^r', with ^ = ^ ^ ^ (B2) 

Zl + Z2 
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and 



4 = ;^^; z' = \{zl + zl). (B3) 

Z-^ ~r 

If ^1 = ^2 = -2 then ^ — \ (this is the case for the same initial and final Gaussian wave 
functions) the formula simplifies to: 

g-iu.r'g-i.V'Y dte^"-*e-2^'*'. (B4) 

1. Diagram 1 

In this case the integral to be calculated is 



j dadbidbs^^ (^2a - <P*c (^2a + - 2k^ <Pb (bi, ba, -k - bi - ba) 

X(/.^(bi + k-k',b2,-k-bi-b2). (B5) 

Here 4>a (with Gaussian parameter and 0c (with Gaussian parameter P2) represent the 
two kaon wave functions and 0b (with Gaussian parameter cti) and (po (with Gaussian 
parameter represent the two nucleon wave functions. 

One sees immediately that the kaon wave function integral on a factors off so that we 
could do the kaon integral directly in momentum space but to be uniform in procedure we 
transform to r-space so we have 

/ dae ^^'i oc / ciacir5dr6e<'^-^)--«e-H"+T+^-^)-^«e-5/3?'-ie-^'^2'''i (B6) 

■k-k' 52 2 (k-k')2 

OC e'^'^^e-^ '■s oc e ^^^d+p)^ . (B7) 
For nucleon B the three momenta bi, b2, and — k— bi — b2 give for the Jacoby coordinates 

p = bi-b2, Q = ^i^, q=-k-^(bi + b2) (B8) 

and hence the B nucleon wave function factor becomes 

(bi-b2)2 [k+|(bi+b2)]' 
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For the D nucleoli the coordinates bi + k — k', b2, and — k — bi — b2 give for Jacoby 
coordinates 

p = k-k' + h,-h2,Q= ^'^'\^'^^\ and q=-^k+lk'-^(bi + b2) (BIO) 



and the wave function becomes 



(k-k'+bi-b2)2 [!k-^k'+|(bi+b2)] 



The integrals on bi and b2 give the delta functions 

6{ti - rs + ^rg - ^r4), 6{-ri + rs + ^rg - ^u). 

We use them to eliminate r2 and to find r2 = ri and r4 = ra. 
The integral becomes 



2 3 OC e 8a2 e 2432 q<- g 6a2 



Combining with the result for the kaon we have an exponent of 



(k-k 



/\2 



+ 



n,2 



[4/32(1 + p)2 6a 



A;2(l-/i) 



+ 



[2/32(1 + p)2 3a2J 



(Bii; 



Writing these factors as r-space transforms we have 

j cibicib2cirirfr2rfr3rfr4e*^^^-^^)"ie-*('^-'^'+^^-''^)-^^e-"'('''+''2) (B12) 

^gi[k+|(bi+b2)]-r3g-*[fk-|k'+f(bi+b2)]-r4g-|a2(r.2+r|)_ 



(B13) 



(B14) 



(B15) 



which agrees (in the limit of equal a's and /?'s) with Eq. 41 of Ref. A numerical check 
was performed by doing the integral on the momentum a and then doing the remaining 
six-dimensional integral by quadrature. 

2. Diagram 2 



The integral corresponding to diagram 2 is 

^ dhidcddi(j)A (2c + - 2k' ) 0c I 2c + '^^^ 



1 + P 



1 + P 



(bi,c, -k - bi - c) 
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X 01)(di, k - k' + bi + c - di, -k - bi - c). (B16) 

Note that there is a typographical error in equation 35 of Ref.^^ where there is a minus 
sign in front of the term 2k/ (1 + p) in the kaon A wave function [16]. This is only a error of 
transcription; we agree with all of their final results. 

For nucleon B the Jacoby variables are (we leave out Q which is not needed) p = bi — 
c, q = — k — |(bi + c) and for nucleon D they are p = 2di — k + k' — bi — c, q = 
— |k + |k' — |(bi + c). With the vectors and being associated with the kaon wave 
functions the imaginary part of the exponent arising from the Fourier transform is 

(bi-c).r,-|k+|(bi + c)].r2-|2di-k + k'-b,-cl.r3+||k-ik' + |(bi + c)].r4. (B17) 

The integrals on di, bi and c give = 0, Y2 = |ri + r4, and rg = rs — ri so that the 
integral becomes 

(B18) 

Defining 7f = (3f/{l3f + we can rewrite the rs integral using Eq. IBll as 
Now we do the integral on using Eq. IBlt 

(B20) 

The contribution of the s and t integrals to the exponent is 

4(k-ko^ 1 (k-kp^ 2 _ p /ifc' e i^e 

(l + p)2 I6/32 ^ 4 1252 ~ 2(1 + pyp^ ~ 2(1 + p)2/32 + 1252 ~ 12^ ^ '^^^ 
so the contributions of the s and t integrals to and 62 are 

' ^ + t4; .\... + t4- (B22) 



2(l + p)2/32 1252' 2(l + p)2/52 12a 
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The coefficient of ri is now 

+ 2^^^^ + -7r(k' - k) (B23) 

^ [(2 + 7f )(1 + P) - 672^]k + [-7i"(l + p) + 6(7f + pW ^ ,^ 

3(1 + P) 

wfiere 

7f = (B24) 



Tlie Gaussian in tlie ri integral becomes 



Q-(o'i+^Pt2+3'^i2M (B25) 



so tfiat tfie contribution of tfie ri integral to the exponential is 



2 I 1^,2 \ ^ ^ 2 I 0/52 I 1^2 \- (B26) 



so that the contributions to 02 and &2 are 

/^2 I /^/2 

+ <-^9 1 '-^2<-^2 



A(al + 2f3f, + lal,y ^"""^ 2^ + 2/3?^ + lO ' ^^^^^ 
A numerical check was done by first calculating the integral on d in the "p" factor of 
nucleon D. It can be integrated on freely and gives only a contribution to the normalization 
of 7riQ;2- The rest of the integral was done by quadrature. 

3. Diagram 3 

The integral for diagram 3 is 

J dadh2dc(l)A (2a - "^c + YTp) '^■s ~ + b2, -a - b2 - k') 

X (/);^(a,b2,-a-b2-k'). (B28) 

The Jacoby variables for nucleon B are: p = — k + k' + a — b2; q = |k — |k' — |(a + b2) 
and for nucleon D: p = a — b2; q = — k' — |(a + b2), giving the integrand 

^i|(a-b2-k+kO-ri+i[k-3k'-3(a+b2)]-r2-(a-b2)]-r3+[k'+|(a+b2)]-r4+(2a-|£|)-r5-(2c+ ^^29) 
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_n.2r.2 _3„2 2 _„2 2 _3 2 2 _2rt2 2 _ofl2 2 



Integration on the momenta gives Tq = 0, r3 = ri + r5 and r2 = r4 + |r5. 
With z = k' — k we have 

(B30) 

^ j dridr4dr5e^^"-''i-5-'-4+[^(k-3kO-|f|]-r5}g-a?rfg-ai(ri+rB)2g-|air|g-faf(r4+|r5)2g-2^ 

(B31) 

^ Jdr^J cite*" *e-2"'*'e-*^^"""^e-"?^'^5y" dse-'^^''^!'^'" e^-'^'^-'^'e-^"^'^ (B32) 



xe 



i[i(k-3k')-|f|]TBg-2^M_ 



The s and t integrals give a contribution to the exponent of 



''+4, = 4 = 4-^- (B33) 



8(52 24(52 6(52 3g^2 3^2 

With 



the coefficient of r5 is 

47f - 1 k - 3k' 2pk 



Cgk + C^k'. (B35) 



3 3 1 + p 

The r5 integral is now 

J (ir5e'(^3'"+^3''')-'"5e-(5"?2+2/3?)'-i (B36) 
so its contribution to the exponent is 

4(|a2, + 2/32) 8(2^2, + 3/32) I ) 

and the contributions to as and 63 are 

8{2ah + 3P!) 4(2c.2, + 3/?2)- ^^"'^ 

A numerical check was done by first doing the integral on c which only appears in kaon 
C. It can be integrated on freely and gives a contribution to the normalization of (27r)^/3|. 
The remainder of the integral was done by quadrature. 
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4. Diagram 4 

The integral for diagram 4 is 



J deidh,dc<PA (^2a - <P*c (^2c + 0s (bi, c, -k - bi - c) 

X 0^(k-k'-a + bi + c,a,k-bi -c). (B39) 

For the Jacoby variables for nucleon B we have: p = bi — c; q = — k — |(bi + c) and 
for nucleon D: p = k - k' - 2a + bi + c; q = -|k + ^k' - |(bi + c). 



(B40) 



2 '2 2' 

The imaginary part of the exponent is 

i I (bi - c) • ri - [k + ^(bi + c)] • r2 - (k - k' - 2a + bi + c) • ra 

-l-5u.iu'-^(b..c),...(2a-i^)..-(.c-M.).. 

The a, bi and c integrals give = —r^, Tq = — ri and r2 = |ri + |r5 + r^. 
The Gaussian factors are 

and the imaginary part of the exponent 
-2 Qk+ .r, + i(k-k')-r4+ Qk - k' - Ts = Urri + U4T4 + U5T5. (B42) 

We make the transformation 
s^ri + rs; t = ri -rg; ri ^ {vH + w^s)/{v'^ + w^); r5 ^ v'^{s - t)/{v'^ + w^) {BA3) 

V 

where 

v^^al + 2Pl, w''^al + 20i. (B44) 

Now 

V'^ -\-w 

and the integral becomes 



v'^rl + w^l = —^{vH'' + w^s") (B45) 
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The integral on t can now be done to give a contribution to the exponent of 

V\ui - {V^ + W^) (k' - k)2(l - p)2 



The r4 integral becomes 



/ 



and the v integral gives 
The s integral becomes 



{k'-kf 2 (k'-k)2 



/ 



dse 



12(52 



2 2 
-w u-i +u uc 2 q: 



24(52 



(B47) 



(B48) 



(B49) 



I" I 1 „2 \ „2 



(B50) 



wV + v^Us 2 „ 

_|_ j 



2^2 - ^;2 + g£ + (^;2 + ^2)^a 



k + 



k' 



3(^2 + w'^) 



(B51) 



= Qk + C^k'. 

so the contribution to the exponent from the s integral is 

(C'4k + C^k02 ^ A;2(C| + Cf + 2C4C^/i) 
4 ( JS. + K) 4 ( + laL) 

and the contribution to and 64 are 

C4C4 



CI + C'i 



4(4xS 



and 



2 ( JS- + i«?2) 



(B52) 



(B53) 



1 2 \ 
+ 3"l2j 

For a numerical check, the a integral which only affects kaon A and the "p" factor of 



nucleon D was done. With x = a — we have 



/ 



('^-■fe) (k-k'-2a+b+c)^ 



dae ^"1 e 



(B54) 



3 r pk _i 



^_^(k-k+b+c) 



4 V«i + 2/3?, 

With this first integral done, the remaining 6-dimensional integral was done numerically. 
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